NAVAL  POSTGRADUATE  SCHOOL 

Monterey ,  California 


THESIS 

C75M15 


COAST  GUARD  DRUG  INTERDICTION:   A 

RENEWAL -REWARD  APPROACH  TO  DETERMINE 

OPTIMUM  INVESTIGATION  TIME 

by 

Eric  A.  Copeland 

March  1988 

Thesis  Advisor:      Donald  P.  Gaver 

Approved  for  public  release;  distribution  is  unlimited 


T238781 


CURITY  CLASSIFICATION  OF  THIS   PAGE 


REPORT  DOCUMENTATION  PAGE 


I    REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


lb    RESTRICTIVE   MARKINGS 


i.  SECURITY  CLASSIFICATION  AUTHORITY 


3    DECLASSIFICATION /DOWNGRADING   SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF   REPORT 

Approved  for  public  release; 
Distribution  is  unlimited 


PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5    MONITORING  ORGANIZATION   REPORT  NUMBER(S) 


I.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  Schoo 


6b    OFFICE  SYMBOL 
(If  applicable) 


7a    NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


:.  ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  California  93943-5000 


7b    ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Monterey,  California  93943-5000 


i.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b    OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


ADDRESS  (Gty,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK  UNIT 
ACCESSION   NO 


rtlTLE  (Include  Security  Cssriication)    C0AST    GUARD    DRUG    INTERDICTION: 
APPROACH    TO    DETERMINE    OPTIMUM    INVESTIGATION    TIME 


A  RENEWAL-REWARD 


PERSONAL  AUTHOR(S) 


Copeland,  Eric  A. 


la.  TYPE  OF  REPORT 

Master's    Thesis 


13b    T  ME  COVERED 
FROM  TO 


14 


DATE  OF  REPORT   ( Year,  Month,  Day) 

198  8    MARCH 


15 


PAGE  COUNT 
"3  8 


;  supplementary  notation  The  views  expressed  in  this  thesis  are  those  or  the  author 
and  do  not  reflect  the  official  policy  or  position  of  the  Department  of 
Defense  or  the  U.S.  Government. 


COSATI  CODES 


FIELD 


GROUP 


SUBGROUP 


18   SJJBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Coastguard,  Drug  Interdiction,  Reward- 
Renewal  Model,  Simulation 


>   ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  renewal -reward  model  is  developed  to  predict  the  optimum  amount  of 
time  that  Coast  Guard  personnel  should  spend  investigating  a  vessel  for 
Lllicit  substances.  The  optimal  investigation  time  is  determined  with 
respect  to  three  criteria;  maximizing  the  number  of  arrests,  maximizing 
:he  quantity  of  drugs  confiscated,  and  minimizing  the  quantity  of  drugs 
that  escape  detection.  A  simulation  study  indicates  that  the  optimal 
Investigation  time  is  very  sensitive  to  underlying  distributional 
assumptions.  The  basic  service  system  model  may  have  wider 
application,  i.e.,  to  combat  modelling  where  it  may  be  desirable  to 
Investigate  a  potential  target  to  estimate  its  value  before  committing 
Limited  resources.  An  adaption  of  the  model  may  also  be  of  help  in 
allocating  resources  for  mineral  exploration. 


0    DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
£D  UNCLASSIFIED/UNLIMITED       □  SAME  AS  RPT  □  DTIC  USERS 


21    ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

Prof.    Donald   P.    Gaver 


22b  TELEPHONE  (Include  Area  Code) 

(408)     646-2605 


-:2c 


OFFICE  SYMBOL 

5  5Gv 


DFORM  1473,  84  mar 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

0  U.S.  Goxrnmtnt  Printlnf  0"lct:  1»i»— COS-24. 


Approved  for  public  release;  distribution  is  unlimited, 


Coast  Guard  Drug  Interdiction:   A 

Renewal -Reward  Approach  To  Determine 

Optimum  Investigation  Time 


by 


Eric  A.  ,Copeland 

Lieutenant,  United  States  Navy 

B.S.,  Texas  A&M  University,  1980 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  1988 


ABSTRACT 

A  renewal-reward  model  is  developed  to  predict  the 
optimum  amount  of  time  that  Coast  Guard  personnel  should 
spend  investigating  a  vessel  for  illicit  substances.  The 
optimal  investigation  time  is  determined  with  respect  to 
three  criteria;  maximizing  the  number  of  arrests, 
maximizing  the  quantity  of  drugs  confiscated,  and 
minimizing  the  quantity  of  drugs  that  escape  detection.  A 
simulation  study  indicates  that  the  optimal  investigation 
time  is  very  sensitive  to  underlying  distributional 
assumptions.  The  basic  service  system  model  may  have  wider 
application,  i.e.,  to  combat  modelling,  where  it  may  be 
desirable  to  investigate  a  potential  target  to  estimate  its 
value  before  committing  limited  resources.  An  adaption  of 
the  model  may  also  be  of  help  in  allocating  resources  for 
mineral  exploration. 


111 


err  st  & 


TABLE  OF  CONTENTS 


I.  INTRODUCTION   1 

II.  RENEWAL-REWARD   2 

A.  EXPECTED  DURATION  OF  PATROL  CYCLE:  E[C]   .  .  3 

B.  REWARD  CRITERIA  1 5 

C.  REWARD  CRITERIA  2 7 

D.  REWARD  CRITERIA  3 8 

III.  RESULTS 11 

IV.  SENSITIVITY  ANALYSIS   16 

V.  CONCLUSIONS 21 

APPENDIX  A 23 

APPENDIX  B 30 

REFERENCE  LIST 32 

INITIAL  DISTRIBUTION  LIST  62 


IV 


ACKNOWLEDGEMENTS 

I  would  like  to  thank  Professor  Gaver  and  Professor 
Jacobs  for  their  encouragement  and  assistance  in  completing 
this  thesis. 


I.   INTRODUCTION 

One  duty  of  the  United  States  Coast  Guard  is  the 
interdiction  of  drug  smugglers.  When  Coast  Guard  personnel 
board  a  suspect  vessel  to  search  for  contraband,  a  decision 
must  be  made  regarding  the  amount  of  time  allocated  for  the 
investigation.  If  a  careful,  but  lengthy,  search  is 
conducted,  there  is  a  high  probability  of  finding 
contraband  if  it  is  present.  However,  a  lengthy 
investigation  incurs  a  penalty.  While  Coast  Guard 
personnel  are  searching  one  vessel,  other  shipping,  which 
could  have  contraband  onboard,  is  passing  through  the 
patrol  area  without  being  examined.  At  the  other  extreme, 
if  the  suspect  vessel  is  examined  in  a  cursory  manner,  the 
Coast  Guard  Vessel  (CGV)  will  be  able  to  stop  and  search 
more  shipping  but  there  is  a  higher  probability  of  not 
discovering  the  contraband,  if  present,  because  of  the 
short  investigation  time.  The  purpose  of  this  thesis  is  to 
develop  a  renewal-reward  model  to  determine  the  optimum 
investigation  time. 


II.   RENEWAL-REWARD  MODEL 

Consider  the  following  scenario.  All  shipping 
encountered  by  a  patrolling  CGV  is  stopped  and  searched  to 
determine  if  contraband  is  onboard.  After  leaving  port,  a 
CGV  searcher  a  time  Si  until  a  vessel  is  sighted.  A 
detected  vessel  has  a  probability  Pb  of  being  Bad  (having 
contraband  onboard)  and  PB  =  1  -  Pb  of  being  Good  (having 
no  contraband  onboard) .  A  Bad  vessel  contains  the  random 
number  J  =  j (1 , 2 ,  3, . . . )  units  of  illicit  substance  with 
probability  bj  .  The  time  required  to  find  the  first  unit 
of  the  contraband,  given  that  J  units  are  present,  is 
modelled  here  as  the  minimum  of  J  independent  identically 
distributed  times  having  distribution  function  Fz  (t) ;  this 
model  is  illustrative  only,  and  may  be  altered  in  various 
realistic  directions.  If  time  T  is  required  to  locate  an 
incriminating  unit,  then 

P{T>T|J=j}  =  [l-Fz(t)]J  ;  j=l,2,3,... 

so,  upon  removing  the  condition, 


1-Ft  (t)  =  P{T>tl  =  S  [l-Fz(t)]Jbj  ;  j=l,2,3,...     (2.1) 


Adopt  the  following  decision  rule:  Establish  a 
predetermined  investigation  time  L.  When  a  vessel  is 
detected,  it  is  stopped  and  searched.   If   no  contraband  is 
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discovered   in   the   investigation   time   L,   the  vessel  is 

released  and  the  CGV  resumes  patrol.   If   any  contraband  is 

discovered  before   the  end  of  the  investigation  time  L,  the 

vessel  is  detained  and  escorted  to  port  for   further  search 

and  investigation.    A   time  period  D  is  required  to  escort 

the  vessel   to   base.     The   patrol   cycle   ends   when  the 

escorted  vessel  arrives  at  the  base. 

The   long   run   average   reward   per   cycle,   R,  can  be 

calculated  in  the  following  manner  [Ross:  pp.  279-294]. 

E[Rc] 

R  =  (2.2) 

E[C] 

where:   E[Rc]  is  the  expected  reward  for  a  patrol 
cycle. 

E[C]  is  the  expected  duration  of  a  patrol 
cycle . 

The  reward  associated  with  an  apprehension  can  be  defined 
in  various  ways.  Three  different  reward  criteria  will  be 
examined.  For  each  of  the  three  cases,  we  will  determine 
the  optimal  investigation  time  L  which  maximizes  the  long 
run  average  cycle  reward. 

A.   EXPECTED  DURATION  OF  PATROL  CYCLE:  E[C] 

We  must  distinguish  between  the  following  two  cases  for 

the  investigation  time  I: 

Ig  -  an  investigation  time  that  results  in  the  release  of 
the  vessel. 

lb  -  an  investigation  time  that  results  in   the  detention 
of  the  vessel. 


Where 


I   = 


Ifl    =   L        with   probability   Pg    +   Pb Ft  (L)    =   a 
lb    =   T        with   probability   Pb  Ft  (L)    =   a" 


The    length   of    a   CGV  patrol    cycle    can   be   represented   as: 


C  = 


Si  +  lb  +  D   With  probability  <x(L) 
Si  +  Ifl  +  C*   With  probability  <x(L) 


where  C*  has  the  same  unconditional  distribution  as  C,  for 
if  the  first  investigation  results  in  a  release  of  the 
vessel  being  searched,  the  process  re-starts  (regenerates) . 
The  expected  length  of  a  patrol  cycle  can  be  expressed 
as : 


E[C]  =  E[S]  +  (E[Ib]  +  E[D])a  +  (L  +  E[C'])a 


or 


E[c]    =   E[S]    +   E[D]a   +   E[Ib]a   ± 

1    -   a 


La 


(2.3) 


(2.4) 


which   can   be   written   as 


[C]    =   E[D]    +   E[Ib]    +   ^fiM-  +   L= 


(2.5) 


Solving  for  the  term  E[Ibl  yields 

P  tFT  (dt) 
E[Ib]  =  °   Ft  (L) 
Integration  by  parts  shows  that 


(2.6) 


P  Ft  (t)dt  =  Ft(L)L  +  [L  tF(dt) 


(2.7) 


Hence 


E[Ib]  = 


i: 


Ft  (t)dt  -  LFt  (L) 


(2.8) 


Ft  (L) 


Therefore 


E[S]     +    E[D]PbFT(L)     +    Pb  FT(t)dt    +    LP 


i: 


E[C]  = 


PbFT  (L) 


(2.9) 


which  reduces  to 


E[C]    =    E[D]    +      ECS]  + 

PbFT  (L) 


o  FT(t)dt  +    P.L 


(2.10) 


Ft  (L)     PbFT  (L) 


B.   REWARD  CRITERIA  1 

Suppose  the  CGV  is  rewarded  for  making  an  arrest. 
Since  only  one  arrest  is  made  per  cycle, 


E[RC  ]  =  1 


(2.11) 


and  from  equation  2.2 


R  = 


E[C] 


(2.12) 


Consider    the    following   example 

F2  (t)    =    1    -    e"Ut  ,    u   =    0.08 

S(t)    =    1    -    e"Xl t  ,    Xi    =    1.1 

D(t)    =    1    -    e"X2t  ,    X2    =    0.3 

bj     =    (1    -    3)3j_1  ,    B    =    0.8;    j    =    1,2,... 

Pb    =    0.2 

From   equation    2.1 

o,  e-n*  (1-B) 

Ft  (t)    =    -Z.     (e-w*)J  (l-B)BJ-*     =    (2.13) 

j=l  1    -    Be""1 


and 


!    _    e-ut 

Ft  (t)     =    1    -    Ft  (t)     =    (2.14) 

1    -    Be-"* 


From   equation   2.10,   the   expected  duration  of  the  patrol 
cycle  is 


1 
E[C]  =  —  + 


•A,  (wAXw*"l/w) 


X2     Pt    1-e-vt 


A-Be-"1-)     f,1   e"UL/l-3e-wuJ 

PflL 

(2.15) 


(1-*-/i-P^.l) 


Pb 


and  from  equation  2.12,  the   long   run   average   reward  per 
cycle  can  be  calculated. 


lXiX2Pb  (2.16) 

R  =  \1-Pe-* 

1-e-wu  \  /l-P\  /l-Pe-«"- 

Pb  (    Xi    +  X2    +   XiXzPb      Ln        )+   P9LXi\2 

1-Pe-uLy  \    Pp/  \       1_P 


C.   REWARD  CRITERIA  2 

The  reward   assigned  to   the  CGV   is  the   amount  of  drugs 
confiscated. 

E[Rc ]  =  E[J|T  £  L]  (2.17) 

The  probability  that  a  vessel  being  searched  contains  J  = 
j  units  of  drugs  and  that  the  drugs  are  discovered  is 

P{(T<JLr(J=j)  }  =  ["l  -  (1-Fz  (L))J~|bj     (2.18) 


so 


P{T£L|«Ft(L)  =  2Z  [l  -  (1-Fz  (L)  )J~|bj      (2.19) 

3=1  L  J 


and  thus 

(1-Fz  (L)  )J~|bj 


['- 


and 


PHJ-j)  I  (T£L)  )  =  (2.20) 

Ft  (L) 


S  j[l  -  (1-Fz  (L)  H~|bj        (2.21) 
j=l   L  J 

E[J|T«SL]  =  

Ft  (L) 


using  the  parameters  from  the  Example  in  Model  1, 

S  j[l  -  (e-»L  )j]  (l-p)pj-1 

j-1    L  J 

E[J|F£L]  =  (2.22) 

(i -e-uL )  / 

*x  e    V (l-3e-»L ) 

which  can  be  evaluated  as  in  the  previous  example  as 

E[J|T£L]     =    

(1-P»[(V   (1-P))     -     (""l /«!-»•■»))] 

(2.23) 

using   equation   2.2,    the    long   run   average   reward   is 

_   E[J   T*LJ 
K   "         E[C] 

or  (2.24) 


R   = 


[(— H-^— )] 

L  \(1-P)V       \  (l-pe-"L  )V   J 
L\(l-P)/       \(l-pe-»«-)/     J 


i        f— vn/ 


1  Xi  \    &u/        \         1-3    /  PcL 

—     +  —   +  +     


X2  B    /  l-e-»"-    \  /  l-e-fL      \  n     /l-e-"L 

Pb   /    \  /     \  Pb    /    


\l-pe-wL/  \  l-ee-"L   /  Xl-pe-wL/ 

D.   REWARD  CRITERIA  3 

The  reward   examined  in  model  three   is  the  difference 
between  the  quantity  of  drugs  confiscated   and  the  quantity 
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of  drugs   that  escape   detection.   The  expected  reward  in  a 
cycle  is 

Re  =  E[J|T£L]  -  E[M] 

where  M   is  defined   as  the   quantity  of   drugs  that  escape 
detection  in  a  cycle. 


M   =    < 


Jib    +    Jd  With   probability    Pb Ft  (L) 

Jig    +    Jo    +   M         With   probability    Pfl     +    Pb Ft  (L) 


where: 


Jib  =  Drug  quantity  passing  through  patrol  area  while  CGV 
searches  a  Bad  ship 

Jd   =  Drug  quantity  passing  through  patrol  area  while  CGV 
escorts  a  Bad  ship  to  shore 

Ji  b  =  Drug  quantity  passing  through  patrol  area  while  CGV 
searches  a  Good  ship 

Jo    =  Drug  quantity  onboard   a  ship   classified  as  Good 
that  is  actually  Bad 


If  conditional  expectations  are  taken  assuming  that  the  Bad 
vessels  pass  through  the  region  according  to  a  Poisson 
process  with  rate  Xi Pb ,  then 

E[M]  »  XiPbECJ]  (E[Ib]+E[D]  )PbFT  (L)  +  (XiPbE[J]E[I9  ]+E[M]  ) 
(Pfl+PbF7(L)  )+E[Jo  ]PbFT(L)  (2.25) 


or 


E[M]    =. 


XiPbE[J]   <   Pb  Ft  (L)dt-LFL 


i: 


J+E[D 


]PbFT  (L)+L       Pg+PbFr  (L 


PbFt  (L) 


E[Jo]PbFT  (L) 


PbFi  (L) 


where   E[Jo]    is   defined   as 


>} 


(2.26) 


j(l-Fz  (L))Jbj 


E[J|T>L]     = 


j=l 


Ft  (L) 


Using  the  distributions  defined  in  the  Example  for  Model  1 


PbX] 


E[M]  = 


r  1-p    /l-pe-»L\     1-e-*11- 
Lnf  ]+  


(1-P) 


1-e-uL 


LP 


o  -t 


3u 


1-3  /  A2  (1-pe-wL )    Pb 


l-^e-i"-  / 
(l-p)e-uL 


(l-3e-ML)  (i-e-ui) 


(2.27) 


and  the  long  run  average  reward  can  be  expressed  as 


R  =  Equation  [2.21]  -  Equation  [2.27]    (2.28) 
Equation  [2.15] 
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III.   RESULTS 

The  optimum  investigation  time  was  determined  for  the 
above  example  for  each  reward  criteria  using  two  methods. 
The  first  method  consisted  of  writing  a  Monte  Carlo 
simulation  of  the  patrol  cycle  using  the  distributions  of 
the  example  and  running  the  simulation  using  various  values 
of  the  investigation  time  L.  Statistics  were  gathered 
during  the  simulation  allowing  the  calculation  of  the  long 
run  average  reward  R.  The  investigation  time  was  varied 
from  zero  to  six  hours  in  increments  of  two  tenths  of  an 
hour.  The  optimum  investigation  time  was  determined  by 
graphing  the  long  run  average  reward  as  a  function  of 
investigation  time  and  finding  the  value  of  L  which 
maximized  R.  Each  cycle  was  replicated  20,000  times.  A 
detailed  discussion  of  the  simulation  can  be  found  in 
Appendix  A. 

THe  second  method  consisted  of  writing  a  computer 
program  for  the  three  equations  representing  the  long  run 
average  reward  for  the  three  reward  criteria  for  the 
example  and  solving  the  equations  for  various  values  of  the 
investigation  time  L.  As  in  the  case  of  the  simulation, 
the  value  of  L  was  varied  from  zero  to  six  hours  in 
increments  of  two  tenths  of  an  hour.  The  optimum 
investigation  time  was   again   determined   by   graphing  the 
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long  run  average  reward  as  a  function  of  investigation  time 
and  finding  the  value  of  L  which  maximized  R. 

Figures  1,2,  and  3  contain  the  results  of  both  the 
simulation  and  numerical  solution  using  the  distributions 
of  the  examples  for  reward  criteria  1,2,  and  3 
respectively.  The  solid  line  represents  the  analytical 
solution  and  the  circles  are  the  simulation  results.  The 
maximum  long  run  average  reward  using  criteria  1  can  be 
achieved  by  using  an  investigation  time  between  1.4  and  2.2 
hours.  The  maximum  long  run  average  reward  using  criteria 
2  and  3  can  be  achieved  using  an  investigation  time  between 
1.2  and  1.4  hours.  It  is  interesting  to  note  that,  using 
the  distributions  presented  in  these  examples,  and 
investigation  time  exists  that  maximizes  all  three  reward 
criteria  simultaneously.  This  occurs  at  1.4  hours  for  the 
input  distributions.  We  do  not  know  that  this  state  of 
affairs  will  persist. 
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IV.   SENSITIVITY  ANALYSIS 

In  order  to  check  the  robustness  of  the  model  with 
respect  to  the  underlying  assumptions,  the  simulation  was 
run  exactly  as  before  with  the  exception  that  Fz  (t)  is  now 
assumed  to  have  a  lognormal  distribution  instead  of  an 
exponential  distribution.  Each  cycle  was  replicated  20,000 
times.  Three  separate  cases  are  examined.  In  case  A,  the 
distribution  function  Fz  (t)  has  the  same  mean  and  variance 
as  the  exponential  distribution  used  in  the  previous 
examples.  In  case  B,  the  distribution  function  Fz  (t)  has 
the  same  mean  but  twice  the  variance  as  the  exponential 
distribution  used  in  the  previous  examples.  Finally,  in 
case  C,  the  distribution  function  Fz  (t)  has  the  same  mean 
but  four  times  the  variance  as  the  exponential  distribution 
used  in  the  previous  examples. 

Figures  4,5,  and  6  contain  a  comparison  between  the 
exponential  base  case  and  the  lognormal  cases  mentioned 
above  for  reward  criteria  1,2,  and  3  respectively.  It  is 
readily  apparent  that  the  results  for  the  exponential  base 
case  and  lognormal  case  A,  with  the  same  mean  and  variance, 
are  different.  Furthermore,  within  the  lognormal  family  of 
curves,  it  can  be  seen  by  examining  case  B  and  C  that  the 
results  vary  significantly  as  the  variance  increases. 
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In  order  to  better  understand  this  behavior,  it  is 
helpful  to  compare  the  quantiles  of  T,  the  time  required  to 
locate  an  incriminating  unit  given  that  J  units  are 
present,  for  the  four  distributions. 


Quantile 

Base  Case 

Case  A 

Case  B 

Case  C 

0.10 

0.242 

1.604 

0.908 

0.431 

0.25 

0.852 

2.502 

1.467 

1.681 

0.50 

2.465 

3.945 

2.856 

1.681 

0.75 

6.504 

6.858 

5.360 

3.909 

0.90 

13.488 

12.155 

10.333 

11.245 

It  can  be  seen  that  some  distributions,  in  particular 
the  lognormal ,  cases  a,  b,  and  c,  have  a  far  greater 
likelihood  of  producing  large  values  of  T  relative  to  the 
others,  such  as  the  exponential  distribution  used  in  the 
examples.  Since  the  time  required  to  find  the  first  unit 
of  drugs  greatly  influences  the  cycle  length,  distributions 
generating  larger  values  of  T  will  produce  significantly 
different  results.  To  further  verify  the  accuracy  of  the 
distributions  produced  by  the  simulation  and  the  quantiles 
listed  above,  Appendix  B  contains  the  calculated  quantiles 
for  lognormal  case  A  produced  using  analytical  distribution 
theory. 
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V.   CONCLUSIONS 

It  has  been  demonstrated  that  a  renewal-reward  approach 
to  modelling  the  Coast  Guard  drug  interdiction  process  is 
feasible  and  that  it  is  possible  to  determine  an  optimal 
investigation  time.  By  considering  several  different 
reward  criteria  in  the  model,  it  is  possible  to  quantify 
the  costs  with  respect  to  the  rest  of  the  reward  criteria 
when  one  criteria  is  selected  and  used  to  arrive  at  an 
optimum  investigation  time. 

Of  the  two  methods  presented  to  obtain  the  optimum 
investigation  time,  the  simulation  approach  is  the  most 
flexible.  The  example  distributions  and  assumptions  used 
in  this  thesis  were  selected  because  they  allowed  an 
analytical  solution  to  be  found.  This  allowed  a  comparison 
of  results  between  the  simulation  and  numerical  solutions, 
thus  verifying  the  simulation.  As  the  model  assumptions 
are  changed  and  different  probability  distributions 
incorporated,  the  simulation  can  be  easily  modified  to 
reflect  these  changes  whereas  an  analytical  solution  may  no 
longer  be  possible. 

As  the  model  is  currently  formulated,  its  usefulness  is 
questionable  due  to  the  sensitivity  to  the  underlying 
distribution  of  the  time  required  to  find  drugs  onboard  a 
vessel  given  drugs  are  present.     As   demonstrated   in  the 
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sensitivity  analysis  even  within  a  given  family  of 
distributions,  the  model  is  sensitive  to  the  distribution 
parameters.  In  order  for  the  model  to  produce  realistic 
results,  rigorous  data  analysis  must  be  conducted  to 
properly  identify  this  distribution. 

The  basic  service  system  model  may  have  wider 
application,  e.g.,  to  combat  modelling,  where  it  may  be 
desirable  to  investigate  a  potential  target  to  estimate  its 
value  before  committing  limited  resources.  An  adaption  of 
the  model  may  also  be  of  help  in  allocating  resources  for 
mineral  exploration. 
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APPENDIX  A 

The  simulation  is  written  in  Fortran  77.  It  consists 
of  a  main  program  and  seven  subroutines.  Uniform  [0,1] 
random  numbers  are  provided  by  calling  procedure  GGUBFS 
located  in  the  IBM  IMSL  single  precision  library.  All  real 
variables  are  computed  using  double  precision  to  minimize 
rounding  error.  Output  is  directed  to  three  separate 
units.  Unit  2  contains  the  calculated  rewards  for  each 
increment  of  investigation  time  using  reward  criteria  1,  2, 
and  3.  Unit  3  contains  the  cycle  length,  drug  quantity 
confiscated,  and  quantity  of  drugs  missed  for  each 
increment  of  investigation  time.  Unit  4  contains  detailed 
information  regarding  the  cycle  for  each  increment  of 
investigation  time.  This  information  includes  the  total 
number  of  ships  searched,  the  number  of  ships  that  are 
good,  bad,  bad  but  declared  good,  and  bad  identified  as 
bad.  Also  included  in  unit  4  output  is  the  time  required 
to  discover  the  drugs  on  a  ship  declared  bad,  the  quantity 
of  drugs  onboard  ships  passing  through  the  area  while  the 
CGV  is  searching  a  ship,  the  quantity  of  drugs  missed  due 
to  short  investigation  time,  and  the  quantity  of  drugs 
confiscated  on  a  bad  ship. 

The   main   program   controls   the   starting   and   final 
investigation  time,  investigation  time  increment  width,  and 
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the  number  of  cycles  per  investigation  time  to  be 
simulated.  Cycle  averages  are  computed  and  output  directed 
to  the  three  units  discussed  above. 

Subroutine  CYCLE  simulates  one  patrol  cycle  and  records 
all  the  relevant  statistics  during  the  cycle.  Subroutine 
STIME  generates  the  random  search  times  required  to  find  a 
ship  from  an  exponential  distribution.  Subroutine  DTIME 
generates  the  random  times  required  to  escort  a  ship  back 
to  base  from  an  exponential  distribution.  Subroutine 
DRGQTY  generates  the  random  quantity  of  drugs  on  a  BAD  ship 
from  a  geometric  distribution.  Subroutine  CLASS  determines 
the  classification  of  a  vessel  based  on  the  deterministic 
value  Pb .  Subroutine  OPCOST  determines  the  quantity  of 
drugs  missed  onboard  other  ships  passing  through  the  area 
while  the  CGV  is  investigating  the  current  ship. 
Subroutine  RTIME  generates  the  random  times  required  to 
find  the  first  unit  of  drugs  given  that  j  units  are 
present. 

Four  versions  of  this  subroutine  are  listed;  one  for 
the  exponential  distribution  and  the  other  three  for  the 
lognormal  case  A,  B,  and  C  distributions. 
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FILE:  DPSIM    FORTRAN   Al 

C     THIS  PROGRAM  SIMULATES  PATROL  CYCLES  FOR  VARIOUS  LENGTHS  OF 
C     INVESTIGATION  TIMES  AND  CALCULATES  THE  VALUES  FOR  THE  THREE 
C     REWARD  CRITERIA.   FUNCTION  GGUBFS  FROM  THE  IMSLSP  LIBRARY  IS 
C     USED  TO  GENERATE  UNIFORM  (0,1)  NUMBERS. 
C 

INTEGER  I,K,MAXT,NCYCLE,Z 

REALX8  L,C,J,M,INC,EC(500),EJ(50Q),EM(500),ITIME(500),TOEND, 
S  R1(500),R2(500),R3(500),TS,GS,BS,BSG,BSB,EIB,OP,MISS,EIG, 
a  CATCH, ES, ED, BJ,DSEED 

DSEED  =  995317. DO 

MAXT  =  <\ 

NCYCLE  =  100 

INC  =  0.20 

L  =  0 

TOEND  =  MAXT/INC 

DO 


200 

K    =    1, TOEND 

L    = 

L    +    INC 

C   = 

0 

J    = 

0 

M    = 

0 

PRINTS, 'COMPUTING  L  =  ',L 
DO  100  I  =  1, NCYCLE 

CALL  CYCLE(DSEED,L,C,J,M,TS,GS,BS,BSG,BSB,EIB,OP,MISS, 
&  CATCH, ES, ED, BJ,EIG) 

100     CONTINUE 

Z=NCYCLE 

WRITER,*)  L 

WRITER, 523)  TS/Z,GS/Z,BS/Z,BSG/Z,BSB/Z 

HRITE(4,524)  EIB/Z, BJ/Z, OP/Z,MISS/Z,CATCH/Z 

WRITE(4,525)  EIG/Z, ES/Z, ED/Z 

WRITE(4,X) 

523  F0RMAT(1X,5F9.4) 

524  F0RMAT(1X,5F9.4) 

525  F0RMAT(1X,3F9.4) 

TS  =0 

BS  =0 

BSG  =  0 

BSB  =  0 

GS  =0 

EIB  =  0 

MISS=0 

CATCH=0 

ES  =0 

ED  =0 

BJ  =0 

OP  =0 

EIG=0 
ITIME(K)  =  L 
EC(K)  =  ONCYCLE 
EJ(K)  =  J/NCYCLE 
EM(K)  =  M/NCYCLE 
RKK)  =  1/ECCK) 
R2(K)  =  EJ(K)/EC(K) 
R3(K)  =  (EJ(K)-EM(K))/EC(K) 
NRITE(2,19)  ITIME(K),R1(K),R2(IC),R3(K) 
WRITE(3,19)  ITIME(K),EC(K),EJ(K),EM(K) 
19       FORMAT  (1X,F6.2,3F9.3) 
200   CONTINUE 

STOP  'PROGRAM  COMPLETE' 
END 

SUBROUTINE  CYCLE(DX,LSTAR,CSTAR, JSTAR,MSTAR,TS,GS, BS, BSG, BSB, 
S  EIB,OP,MISS,CATCH,ES,ED,BJ,EIG). 

REALX8  LSTAR,CSTAR,JSTAR,MSTAR,LAMDA1,LAMDA2,MU,BETA,S,R,QTY, 
a  D, LEAK, TIME, PBAD,TS,GS, BSG, BSB, EIB, OP, MISS, ES, ED, B J, BS, CATCH, 
a  EIG,DX 

INTEGER  NC 

CHARACTERS*  TYPE 

LAMDA1  =1.10 

LAMDA2  =0.3 
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FILE:  DPSIM     FORTRAN   Al 


MU  =  0.08 
BETA  =0.8 
PBAD  =  0.2 

20  TS=TS+1 

CALL  STIMECDX,LAMDA1,S) 
ES=ES+S 

CALL  CLASSCDX, PBAD, TYPE) 
IF(  TYPE  .EQ.  'GOOD')  THEN 
EIG  =  EIG+LSTAR 
GS=GS+1 
QTY=0 

TIME  =  LSTAR 

CSTAR  =  CSTAR  +  TIME  +  S 

CALL  OPCOST(DX,LAMDAl, TIME, BETA, PBAD, LEAK) 
MSTAR  =  MSTAR  +  LEAK 
OP=OP+LEAK 
GOTO  20 
ELSE  IF(  TYPE  . EQ .  'BAD')  THEN 
BS=BS+1 

CALL  DRGQTY(DX,BETA,QTY) 
BJ=BJ+QTY 

CALL  RTIME(DX,MU,QTY,R) 
IF(R  .GT.  LSTAR)  THEN 

EIG=EIG+LSTAR 

BSG=BSG+1 

TIME  =  LSTAR 

CSTAR  =  CSTAR  +  TIME  +  S 

CALL  OPCOST(DX,LAMDAl, TIME, BETA, PBAD, LEAK) 

MSTAR  =  MSTAR  +  LEAK  +  QTY 

OP=OP+LEAK 

MISS=MISS+QTY 

GOTO  20 
ELSEIF(R  .LE.  LSTAR)  THEN 

EIB=EIB+R 

BSB=BSB+1 

CALL  DTIME(DX,LAMDA2,D) 

ED=ED+D 

TIME  =  R  +  D 

CSTAR  =  CSTAR  +  TIME  +  S 

JSTAR  =  JSTAR  +  QTY 

CALL  OPCOST(DX,LAMDAl, TIME, BETA, PBAD, LEAK) 

MSTAR  =  MSTAR  +  LEAK 

OP=OP+LEAK 

CATCH=CATCH+QTY 

GOTO  21 
ELSE 

STOP  'ERROR  1« 
ENDIF 
ELSE 

STOP  'ERROR  2' 
ENDIF 

21  RETURN 
END 


SUBROUTINE  STIME( DRSEED, INPUT, OUTPUT) 

REALX8  INPUT, OUTPUT, RV, DRSEED 

RV=GGUBFS(DRSEED) 

OUTPUT  =  LOG(RV)/(-INPUT) 

RETURN 

END 


SUBROUTINE  DTIMECDRSEED, INPUT, OUTPUT) 

REALX8  INPUT, OUTPUT, RV, DRSEED 

RV  =  GGUBFS( DRSEED) 

OUTPUT  =  LOG(RV)/(-INPUT) 

RETURN 

END 
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FILE:  DPSIM     FORTRAN   Al 


SUBROUTINE  DRGQTY( DRSEED, INPUT, OUTPUT) 

REAL*8  DRSEED, INPUT, OUTPUT, RV 

RV  =  GGUBFS(DRSEED) 

OUTPUT  =  AINT(LOG(RV)/(LOG(INPUT)))+l 

RETURN 

END 


SUBROUTINE  CLASSC DRSEED, INPUT, OUTPUT) 
REAL*8  DRSEED, INPUT, RV 
CHARACTER**  OUTPUT 
RV  =  GGUBFS(DRSEED) 

IF(RV  .LE.  INPUT)  THEN 
OUTPUT  =  'BAD' 

ELSE 

OUTPUT  =  'GOOD' 

ENDIF 
RETURN 
END 


SUBROUTINE  OPCOST( DRSEED, INPUT1 , INPUT2, PI , P2, OUTPUT) 

REAL*8  DRSEED, INPUT1 , INPUT2, PI, P2, OUTPUT, STOR,NSHIP 

OUTPUT  =  0 

JO  =  1/C1-P1) 

NSHIP  =  INPUT2XINPUT1 

OUTPUT  =  (P2*NSHIP)*J0 

RETURN 

END 

SUBROUTINE  RTIME( DRSEED, INPUT1 , INPUT2, OUTPUT) 
:      EXPONENTIAL  DISTRIBUTION  -  BASE  MEAN  8  VARIANCE 
REALX8  DRSEED, INPUT1 , INPUT2, OUTPUT, RV 
RV  =  GGUBFS(DRSEED) 
OUTPUT  =  L0G(RV)/(-INPUTl*INPUT2) 
RETURN 
END 

SUBROUTINE  RTIMEK DRSEED, INPUT1 , INPUT2, OUTPUT) 
;      LOGNORMAL  DISTRIBUTION  -  BASE  MEAN  8  VARIANCE 

REAL*8  INPUT1,INPUT2, OUTPUT, STDV, AVG, PI, LOW, U1,U2,GG, 
8  DRSEED 
INTEGER  II, FINISH 
GG  =  LOGC2.0) 
LOW  =  999999 
FINISH  =  INT(INPUT2) 
STDEV  =  SQRT(GG) 
AVG  =  -(LOG(INPUTl))-  0.5*GG 
PI  =  3.141592654 
DO  662  II  =  1, FINISH 
Ul  =  GGUBFS(DRSEED) 
U2  =  GGUBFS(DRSEED) 
RV  =  SQRT(-2XL0G(U1))XC0S(2*PIXU2)*STDEV+AVG 

RVT  =  EXP(RV) 

IF(RVT  .LT.  LOW)  THEN 

LOW  =  RVT 
ENDIF 
662   CONTINUE 

OUTPUT  =  LOW 

RETURN 

END 

SUBROUTINE  RTIME2(DRSEED, INPUT1 , INPUT2, OUTPUT) 
:      LOGNROMAL  DISTRIBUTION  -  BASE  MEAN  8  2XVARIANCE 

REALX8  INPUT1,INPUT2, OUTPUT, STDV, AVG, PI, LOW, U1,U2,GG, 
8  DRSEED 
INTEGER  II, FINISH 
GG  =  LOG(3.0) 
LOW  =  999999 
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FINISH  =  INTCINPUT2) 

STDEV  =  SQRT(GG) 

AVG  =  -(LOGCINPUT1))-  0.5*GG 

PI  =  3.141592654 

DO  662  II  =  1, FINISH 

Ul  =  GGUBFS(DRSEED) 

U2  =  GGUBFS(DRSEED) 

RV  =  SQRT(-2*L0G(U1))*C0S(2*PI*U2)XSTDEV+AVG 

RVT  =  EXP(RV) 

IF(RVT  .LT.  LOW)  THEN 

LOW  =  RVT 
ENDIF 
662   CONTINUE 

OUTPUT  =  LOW 

RETURN 

END 

SUBROUTINE  RTIME3( DRSEED, INPUT1 , INPUT2, OUTPUT) 
C      LOGNROMAL  DISTRIBUTION  -  BASE  MEAN  S  4XVARIANCE 

REAL*8  INPUT1 , INPUT2, OUTPUT, STDV, AVG, PI , LOW, Ul , U2,GG, 
5  DRSEED 
INTEGER  II, FINISH 
GG  =  L0GC5.0) 
LOW  =  999999 
FINISH  =  INTCINPUT2) 
STDEV  =  SQRT(GG) 
AVG  =  -(LOG(INPUTl))-  0.5XGG 
PI  =  3.141592654 
DO  662  II  =  1, FINISH 
Ul  =  GGUBFS(DRSEED) 
U2  =  GGUBFS(DRSEED) 

RV  =  SQRTC-2*L0G(U1))XC0S(2XPIXU2)XSTDEV+AVG 

RVT  =  EXP(RV) 

IF(RVT  .LT.  LOW)  THEN 

LOW  =  RVT 
ENDIF 
662   CONTINUE 

OUTPUT  =  LOW 

RETURN 

END 
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FILE:  TSTNUM    FORTRAN   Al 

C      THIS  PROGRAM  COMPUTES  THE  VALUES  OF  THE  THREE  REWARD  CRITERIA 
C      FOR  VARIOUS  VALUES  OF  L  FOR  COMPARISON  WITH  SIMULATION  RESULTS 
C 

REAL    PB,LAMDA1,LAMDA2,BETA,MU,INC,ES,ED,FT,FTBAR,A,ABAR, 

5  L,TOEND,EC(100),EJC100),EM(100),R1(100),R2(100),R3(100), 

6  ITIME(100),PG,B,BB,E,C1,C2,C3,C4,M1,M2,TEMP,J1,RATE 
INTEGER   MAXT,K 

MAXT  =  6 
INC  =  0.20 
PB  =  .20 
PG  =  1-PB 
LAMDA1  =1.10 
LAMDA2  =0.3 
ES=1/LAMDA1 
ED=1/LAMDA2 
BETA  =  0.8 
MU  =  0.08 
B=BETA 

DD-1  _  O 

RATE  =  (LAMDA1*PB)/BB 

TOEND  =  MAXT/INC 

L  =  0 

DO  100  K  =  1, TOEND 

L=L+INC 

E=EXP(-MU*L) 

TEMP=1-(B*E) 

FTBAR  =  (BB*E)/TEMP 

FT=1-FTBAR 

ABAR=PBXFT 

A=1-ABAR 

Cl=  ES/ABAR 

C2  =  ED 

C3=(PB/ABAR)X((BB/(BXMU))XL0G(TEMP/BB)-(LXFTBAR)) 

C4=(LXA)/ABAR 

J1=(BB/FT)X((1/(BBXBB))-(E/(TEMPXTEMP))) 

M1=(C3+ED)XRATE+((LXAXRATE)/ABAR) 

M2=BB*(E/(TEMPXTEMP))X(1/ABAR)XPB 

ITIME(K)  =  L 

EC(K)  =  C1+C2+C3+C4 

EJ(K)  =  Jl 

EM(K)  =  M1+M2 

RKK)  =  1/ECCK) 

R2(K)  =  EJ(K)/EC(K) 

R3(K)  =  (EJ(K)-EM(K))/EC(K) 

WRITE(2,132)  ITIME(K),R1(K),R2(K),R3(K) 

WRITE(3,132)  ITIMECK),EC(K),EJ(K),EM(IC) 

WRITE(4,X)  L 

WRITE(4,133)  C1,C2,C3,C4 

WRITEC+,134)    Jl 

WRITEC+,135)    M1,M2 

132  FORMAT    ( IX, F6 . 2,  3F9 . 3) 

133  FORMAT    (1X,4F12.5) 

134  FORMAT    (1X,F12.5) 

135  FORMAT    (1X,2F12.5) 
WRITEC4,X) 

100   CONTINUE 

STOP  'PROGRAM  EXECUTION  COMPLETE' 
END 
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APPENDIX  B 

Recall   that    the   probability   distribution   for   the 
quantity  of  drugs  on  a  vessel  given  that  it  is  bad  is: 


Let 


so 


bj  =  (1  -  3) 3J -1 ;   3  =  0.8,  j  =  1,2,3, 


P{T>tl  ■  Ft  (t)  =  Z:  (l-Fz(t))Jbj 


—     _  (l-p)Fz (t) 
Ft  (t)  -   1_(5Fz  (t) 


and 


p(T<t,  _  r    _  d-P)F2  (t)  _    1-Fz  (t 
niStl   x    l-3FT(t)     l-pPT(t 


now  let 


1-Fz  (tp)  =  p 


l-|3Fz  (tP  ) 
which  can  be  written  as 

and 


Fz  (tP)  =  1  -  Fz  (tP)  =  P\1tKt) 

1-pp 
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Assuming  Fz   (t)  is  lognormal  with  mean  =  2.1792  and  standard 
deviation  =  0.8326  as  in  distribution  Case  A 


therefore 


<D 


In  tP  -    u\=   P(l-P) 
1-3P 


tP  =  EXP 


u  +  azP  <  i  -p  > 

l  -  0  P 


where  ZP  <  i  - p  >  is  taken  from  tables  of  the  Standard  Normal 
i  -  p  p 

distribution. 


Using  the  equation  for  tP ,  the  following  table  of  quantiles 
was  generated: 


0.10 
0.25 
0.50 
0.75 
0.90 


1.65 
2.46 
3.95 
6.78 
11.99 


It  can  be  seen  that   these   quantile   values   agree  closely 
with  those  presented  in  the  sensitivity  analysis. 
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